Metal-insulator transition in one-dimensional lattices with chaotic energy sequences 
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We study electronic transport through a one- dimensional array of sites by using a tight binding 
Hamiltonian, whose site-energies are drawn from an chaotic sequence . The correlation degree 
between these energies is controlled by a parameter regulating the dynamic Lyapunov exponent 
measuring the degree of chaos. We observe the effect of chaotic sequences on the localization length, 
conductance , conductance distribution and wave function, finding evidence of a Metal-Insulator 
Transition (MIT) at a critical degree of chaos. The one dimensional metallic phase is characterized 
by a Gaussian conductance distribution and exhibits a peculiar non-selfaveraging. 
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It is well-known that one-dimensional systems with un- 
correlated disorder behave like insulators because their 
electronic states tend to localize 0. For systems whose 
length is larger than the electronic localization length 
the conductance vanishes exponentially. Nevertheless, 
the Anderson localization)^] does not address the effect 
that correlations have on electronic states in disordered 
lattices. Short range correlations have been shown to 
be responsible for low dimensional delocalization behav- 
ior in direct experimental realizations using semiconduc- 
tor superlattices|3, 0, 0. Such is the behavior of the 
one-dimensional random-dimer model, where transitions 
from localized to extended states are found when the 
Fermi energy assumes values belonging to a discrete set 
of energies |g. 

Extended states have been reported for incommensu- 
rate one dimensional systems which can be viewed as 
long range correlated. This is known to occur within 
the Harper model 7] whose the site energies are given 
by Cn — {w / 2) cos{27muj) where lu is irrational and w 
is a parameter. Below a critical value of w extended 
states are obtained. A more general model was intro- 
duced by Griniasty and FishmanQ including the Harper 
model as a particular case, describes a one dimensional 
model where the energy e„ of site n is generated by the 
relation e„ = cos(7rart''). The resulting sequence is de- 
scribed as pseudorandom depending on the value of v and 
a an irrational number. In this model a mobility edge is 
found whereupon the system exhibits a metal insulator 
transition for a threshold value of v. 

Mobility edges in one dimensional systems have also 
been recently reported in non-periodic Kronig-Penney 
potentials correlated over long distances 0. Here it 
was shown how to construct site energy sequences, that 
lead to predetermined localization-delocalization tran- 
sitions. Interestingly, mobility edges have been found 
experimentally for these systems by microwave transi- 



tion measurements on waveguides with inserted corre- 
lated scatterers.lOl- More recently, non-trivial electronic 
transport properties in low dimensional systems such as 
certain proteins, and DNA chains^, 0| have prompted 
new models involving long range correlated disorder [l^. 

In this work, we study electronic transport through a 
quasi-one dimensional system with site energies obeying 
a chaotic sequence characterized by a dynamic Lyapunov 
exponent A, that gives a measure of the degree of chaos 
and short ranged correlation. The model of chaotic se- 
quences used here is of much interest since it addresses 
in a continuos fashion periodic and chaotic sequences of 
varying dynamic Lyapunov exponent with a single con- 
trol parameter z. The latter parameter can in turn be 
related analytically to the dynamic Lyapunov exponent, 
the rate at which information is lost in the sequence. 
Random sequences are then identified with large values 
of the dynamic Lyapunov exponent were no memory of 
the previous values is preserved [l4j|. while small values of 
the exponent relate to correlated sequences. The manner 
of introducing disorder is reminiscent of that proposed 
in ref. 8] discussed above, where a parameter served 
to measure the degree of pseudo-randomness. Pseudo- 
randomness, in this case, is well defined by a sequence 
of numerical tests in ref. |l5l | in relation to Anderson 
localization. 

This work first addresses the chaotic energy sequence 
recursion and the analytic form of the dynamic Lyapunov 
exponent. The methods for obtaining reliably the wave 
function and the localization length as a function of en- 
ergy follow, along with the computation of the conduc- 
tance of the system from the Green's function. We then 
show that the localization length grows faster than the 
length of the system, within a finite energy band, below 
a critical value of the parameter governing the degree of 
chaos z, indicating delocalized electron behavior. Such a 
conclusion is supported by studying the conductance as 
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a function of system size, its asymptotic scaling proper- 
ties, and the localization of the wavefunction. The con- 
ductance distribution in the extended phase is a narrow 
Gaussian for different realization of the chaotic sequence 
but its mean fluctuates more strongly with system size in 
comparison to the distribution width, displaying non-self 
averaging behavior. 

We shall consider a tight binding Hamiltonian with 
nearest-neighbor coupling in order to describe one elec- 
tron moving through a chain of N sites 

H = J2en\n){n\ + V\n){m\ . (1) 

n {n,m) 

where e„ are the site energies and V the nearest neighbor 
coupling. This chain is attached to semi-infinite leads at 
each end, whose site energies are equal to eo, and the 
hoppings are unity. The site energies in the sample are 
generated by the recursive formula (inspired by the re- 
cursion of Nazareno, Gonzalez and Gosta (NGG)|13j| ) 




FIG. 1: Bifurcation map corresponding to the recursive re- 
lation of Eq. |2| The sequence at the threshold oscillates 
between essentially two values (0 and 1), but gradually grows 
to the full interval ( 0,1) . The sequence is only strictly random 
in the z » 1 limit [l^. 



£«+! = 1 — sin^ [z arcsin Y^e^ , 



(2) 



z being the controlling parameter that regulates the dy- 
namic Lyapunov exponent. By using the transforma- 
tion Un = ^ arcsin we obtain a piecewise linear map 
yn+i — fiVn)- One can then find an analytic expression 
for the dynamic Lyapunov exponent by applying the re- 
lation 



A — lim 
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Therefore, for z > 1 the map is chaoticpU^. 

In Fig. ^ we show a bifurcation map corresponding to 
the formula (|2Jl. We drew energy values after a transient 
of 3000 iterations of the recursion formula, which applies 
to every further computation. Different realizations of 
the chaotic sequence are generated by initiating the re- 
cursion with different initial values. 

We can distinguish different ranges of values for the 
parameter z: periodic region, for z < 1, in which the 
sequence converges to one fix point; and chaotic region for 
z > 1, where e„ varies over a fixed interval. Notice that 
while z moves away from z = 1 (value at which we obtain 
an ordered binary array) the e„-values spread until this 
variable assumes almost all values in the interval (0, 1). 

To compute the wave function for the tight bind- 
ing model, we use the band matrix routines described 
in ref.|l6l| using Gaussian elimination and back substi- 
tution. Diagonalizing by this method one can handle 
large one dimensional systems reliably, in order to as- 
sess the localizing properties of the wave function. On 
the other hand, both the localization length and the con- 
ductance are more conveniently computed by renormal- 
ization procedure that generates the Green's Function 



G(e) = [e + ir] — H]~^ through a decimation method[l7|. 
This method consists in the elimination of sites along the 
chain in an iterative fashion, in such a way that we ob- 
tain an effective Hamiltonian which involves the sites 1 
and n -|- 1 in the chain, whose site energies are 



ei(6) = 6i + A+„^(e), 
e„+i(e) = e„+i + A^^i(e). 



(4) 
(5) 



Here A^^^^ (e) is the energy correction to the first site due 

the presence of the other sites in the chain, and A~_,_-^(e) 
is the energy correction to site n + 1. These corrections 
are given in detail in ref. 0, 0| . 

In order to take into account the leads in the calcu- 
lation procedure, we add a self energy correction to site 
energies ei and e„+i^^ whose value within the band 



|e-eo| < 2|F| is S ^^±i^V^ - where eo is 

the site energy of the semi-infinite chain representing the 
incoming and outgoing leads. Once the Green's function 
is calculated, the localization length ^ is given by 
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(6) 



and the conductance by the two-lead Landauer for- 
mula, expressing the conductance in terms of the Green's 
function 1 1 Si. 



g=^(2r(e)|Gi,Arp2r(e)) 



(7) 



r(e) being the imaginary part of the self-energy correc- 
tion described above. 

In FigI21we plot the localization length as a function 
of the energy (referred to the site energies of the semi- 
infinite leads) for three different values of the parameter z 





FIG. 2: The localization length as function of the energy for 
the shown values of z, eo = 1 and the coupling between sites 
set always to V = —1. The horizontal dashed line indicates 
the size of the system. Two bands of states appear in cor- 
respondance with the ordered alternating sequence close to 
z = l+. 



FIG. 3: Conductance g (units 2e /h) eversus size of the sys- 
tem for different values of z. The curves show a threshold 
value in the vicinity of z = 1.085 above which the conduc- 
tance decays exponentially. Each curve is averaged over 5 
different initializing seeds for the recursion generating chaotic 
sequences. 



in the vicinity of z = 1.1. The dashed hne shows the size 
of the system N demonstrating that > N within a finite 
energy range. In this case two relevant energy ranges 
are observed due to the nature of the initially binary 
energy sequence generated (see Fignj. The localization 
length grows linearly with the size of the system in the 
two energy ranges shown. 

The delocalized nature of states within the energy 
bands observed in Fig|21can also be corroborated by com- 
puting the conductance according to Eq|3 Fig|31 shows 
the conductance 5' as a function of system size for differ- 
ent values of z. Each curve is built by averaging over a 
few different initializing seeds for the chaotic sequence, 
avoiding spurious transient effects as described above. 
Below z — 1.085, the conductance is essentially constant 
up to system sizes of iV = 10^". Above the mentioned 
threshold value, the conductance drops exponentially to- 
wards zero, indicating a clearcut transition. 

As a check that we have reached asymptotic behavior 
we have plotted the conductance 5 as a function of the 
parameter z, increasing size until the curves reached a 
fixed limit. These results are shown in Fig^ where no 
further changes are evident as we increase system sizes 
above 10* sites. The figure suggests a clearcut transi- 
tion between extended states and localized states. The 
inset shows a running average of the main curves so as to 
decrease fluctuations. The computation of the wavefunc- 
tion versus the position along the one dimensional chain, 
as explained above Eq^led to the same conclusion. We 
diagonalized systems of up to 1000 sites for different val- 
ues of the parameter z and energy e — eo = —1.5. The 
wavefunction decays exponentially beyond a certain value 
of the control parameter. For z = 1.5 the wavefunction 
is completely localized below 500 sites. 



To show that the average conductance is a good mea- 
sure of the conductance distribution, the full distribution 
of the conductance is obtained as a function of system 
size and initializing seeds. In FigElwe show that the dis- 
tribution, when sampled for different initializing seeds, is 
a narrow Gaussian (see inset) of width 2% of the mean 
conductance value. On the other hand, as the system 
size is changed the whole distribution fluctuates within a 
broader interval of 30% around the mean conductance of 
FigEl These two different fluctuational amplitudes im- 
ply a non-selfaveraging property of the extended phase of 
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FIG. 4: Conductance (units 2e^ /h) versus z for linear chains 
of the sizes shown. The collapse of the curves demonstrates 
that no further changes operate on the curves and an asymp- 
totic limit has been reached. The inset shows the running av- 
erage of the curves that reduces noise and indicates a clearcut 
metal-insulator transition. 
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FIG. 5: Conductance distribution in the metallic phase as 
a function of systems size. Each distribution is generated by 
changing the seed of the chaotic energy sequence, and they 
are well defined Gaussians (see inset). Nevertheless changes 
in length involve larger fluctuations in the mean values. 



the system rather different from a regular metaUic phase. 
The distribution in the critical region, on the other hand 
is broad and requires a much more extensive study in 
itself. 

As the results show, we obtain a metal-insulator tran- 
sition for a finite band of states in a one dimensional 
system with a chaotic sequence of energies. The metallic 
character is assessed from the behavior of the localiza- 
tion length with system size, the scaling of the conduc- 
tance and the behavior of the wavefunction as a function 
of the chaoticity parameter z. Nevertheless, for the en- 
ergy sequence used, correlations are present but are not 
long ranged, so such feature seems sufficient but not nec- 
essary for delocalization to occur in the present sense. 
One might be tempted to test a model where energies 
around and 1 are pertubed randomly in the vinicity 
of those energies, generating a small band of values (re- 
sembling structure of Fig^ Nevertheless, such a model 
behaves asymptotically localized for arbitrarily small dis- 
order. We have thoroughly tested this conclusion with 
the methods described before. 

Recently, several physical systems have been reported 
to be very sensitive to changes in the nature of disorder 
within the system. Among these systems we can men- 
tion: the two-dimensional Ising model, the problem of 
ballistic deposition, and random walks [itA l2Cl |21| . Re- 
garding random walks, it is possible to show that if there 
is "determinism" in the external noise applied to random 
walks, the resulting walks observe symmetries which are 
absent when the noise is completely random. Anomalous 
behavior ensues in the computation of the mean distance 
versus time. 



In the context of ballistic deposition, that serves as 
a prototype for studies of dynamic scaling phenomena 
in nonequilibrium growth processes, the external influ- 
ences are usually modeled by random noise. D' Souza et 
al have studied ballistic deposition using different deter- 
ministic sources of noise jjOj, and find fluctuations sta- 
tistically inconsistent with the steady-state distribution 
and ergodicity i.e. driving the dynamics of the system 
with different deterministic sources of "noise" results in 
selective sampling of the phase space. Such determinis- 
tic noise has been shown to modify the scaling exponents 
predicted by theory assuming strictly random noise. 

The class of energy sequences addressed here seem to 
open new possibilities, which rather than sequence cor- 
relation, emphasize determinism implicit in chaos. Such 
determinism seems to play a non-trivial role in combina- 
tion with the quantum mechanics of localization. This is 
a new example of physical systems where deterministic 
disorder can produce new phenomena. However, within 
quantum mechanics the phenomenon is more subtle: If 
the Lyapunov exponent that defines the chaos of the en- 
ergy sequences is larger than certain critical value, the 
system behaves as the usual Anderson model. 
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